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Abstract 

We review the results of our recent numerical investigations on the electronic properties of disordered two dimensional 
systems with chiral unitary, chiral orthogonal, and chiral symplectic symmetry. Of particular interest is the behavior 
of the density of states and the logarithmic scaling of the smallest Lyapunov exponents in the vicinity of the chiral 
quantum critical point in the band center at E = 0. The observed peaks or depressions in the density of states, the 
distribution of the critical conductances, and the possible non-universality of the critical exponents for certain chiral 
unitary models are discussed. 
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1. Introduction 

Two-dimensional disordered systems have been at- 
tracting special attention for many years because d — 2 
is the lower critical dimension of the metal-insulator 
transition (MIT) [1]. For lattice systems with orthog- 
onal symmetry (random on-site disorder with time re- 
versal symmetry) all electronic states are localized in 
the limit of infinite system size. However, for weak dis- 
order and energies close to the band center, the local- 
ization length can become very large. On length scales 
smaller than the localization length, the wavefunctions 
exhibit self-similar (fractal) behavior yfl. The presence 
of spin dependent hopping changes the symmetry of the 
model to symplectic, and enables the system to undergo 
a metal-insulator transition at a certain value of the dis- 
order strength J31-01 ■ The critical eigenstates at the MIT 
exhibit multifractal properties Jisj] and the localization 
length was reported to show a parity dependence fl. 
A strong magnetic field turns the symmetry to unitary 
and induces critical states 111 OA . i.e., singular energies 
where the localization length of the multifractal eigen- 
states ll ill 1211 diverges, which are important for the ex- 
planation of dissipative transport Il3l [l4ll in the quan- 
tum Hall effect. 



Two-dimensional (2D) models possessing an addi- 
tional chiral symmetry exhibit various electronic prop- 
erties not observed in the situations mentioned above. 
The chiral symmetry can be found in models defined 
on bi-partite lattices with non-diagonal disorder only 
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II15L 11611 . Despite the disorder, the energy eigenvalues 
appear in pairs, E n and -E„ symmetrically around the 
band center E — (for the definition of chiral 2D mod- 
els, see section Chiral symmetry implies unusual 
properties of the model in the vicinity of the band cen- 
ter E = 0. For most chiral cases, the density of states 
and the localization length are di verging and the band 
center is a quantum critical point jl7ill8ll . At zero tem- 
peratures, an infinite sample is metallic at E — but in- 
sulating for any non-zero energy. The appearance of the 
criticality at E = originates from the chiral symmetry. 
However, the existence of the critical point also depends 
on the boundary conditions. As listed in Table Q] the 
sample exhibits chiral symmetry only for special combi- 
nations of boundary conditions and parity. This bound- 
ary and parity dependence of the sample's length L z and 
width L has no analogy in 'standard' disordered models. 

The special symmetry of the energy spectra may be 
accompanied by a non-analytical behavior of the den- 
sity of states (DOS) at the chiral critical point [18, 23- 
12511 . In 2D chiral unitary models defined on a bricklayer 
112611 . which has the same topology as graphene's hon- 
eycomb lattice, the DOS exhibits a sharp drop near the 
band center going to zero at E — and depends on both 
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Table 1: Two dimensional models with non-diagonal disorder possess 
the chiral symmetry only for special choices of the boundary condi- 
tions and the parity in the number of sites. For a given combination of 
boundary conditions, periodic (P) and Dirichlet (D), the chiral sym- 
metry (Ch) and chiral symmetry with an extra eigenvalue at E = 
(Ch+) is observed 



disorder and system size 12711 . Contrary to this behavior, 
the DOS of the chiral orthogonal system is finite at the 
band center and showing a narrow extra peak in the case 
of square lattice samples (see below). 

Similarly to other critical regimes, systems with chi- 
ral symmetry can be analyzed using the single parame- 
ter scaling theory 01. However, the scaling parameter 
is not the ratio of the system size L to the correlation 
length £(£), but the ratio of the logarithm of these pa- 
rameters instead 112811 
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Also, the energy dep endence of the correlation length is 
logarithmic HEIB, 



In(f(£)/f )~|ln(£o/|£|)l*, 



(2) 



in contrast to the power-law scaling dependence tj(E) ~ 
\E - E c \~ v observed in non-chiral disordered systems. 
Thus, models with chiral symmetry enable the detailed 
analysis of logarithmic scaling, discussed previously 
in 10]. 

Chirality also strongly affects the transport proper- 
ties of the system. The non-Ohmic behavior of chiral 
systems with an odd number of open channels, where 
the mean conductance decreases much more slowly 
with the length of the system (see Fig. [TJ, has been 
predicted theoretically 12211 and confirmed numerically 
in EH. 

In this paper we review some results of our recent 
investigations on electric transport properties of two di- 
mensional chiral systems. In the following section we 
briefly introduce the models studied. In section |4] we 
summarize our findings for the Lyapunov exponents, the 
critical conductance and its probability distribution, and 
show recent calculations for the density of states. Of 
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Figure 1 : The length dependence of the mean conductance (g) for 
quasi-one dimensional chiral systems of length L. and fixed width 
L = 65 a. The energy is E = 0. Due to the chiral symmetry of the 
model, the crossover from the Ohmic l/L--behavior (dotted line) to 
the 1 / VZT dependence (dashed line) is observed, in agreement with 
theoretical predictions f22ll . Two models with chiral unitary (ChU: 
/ = 0.25 h/e) and chiral orthogonal (ChO: W/to = 1) symmetry, de- 
fined in section [2] were considered. 



special interest is the scaling analysis of the diverging 
critical electronic states at E — |24, 3^, 31]. 



2. Models 

In the absence of diagonal disorder the single-band 
tight-binding Hamiltonian defined on the sites « of a 
two-dimensional bricklayer or square lattice with lattice 
constant a reads 



< 1~^ — ^ ^ tr 

(ntn') 



(3) 



where the sum is over nearest neighbors only. The 
random disorder is incorporated in the hopping terms, 
which also determine the symmetry of the problem. 
Square lattice and bricklayer lattice differ only in the 
absence of every other vertical bond in the latter and so 
the coordination number is reduced to three (see Fig.[2j. 

2.1. Unitary symmetry 

To describe a disordered chiral 2D system with bro- 
ken time-reversal symmetry, the hopping terms in the 
(transversal) x direction are chosen to acquire complex 
phases and are defined as 



t x = t e" 



(4) 
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where for a bricklayer lattice the phases 9 x ^ x+a , z — 
#i,z+2a;A+a,z+2a - "x^.i-.z are determined by the total flux 
threading the plaquette at (x, z) 



, ph 
<& x , = - - + i 

q e 



(5) 



Here, p and q are mutual prime integers and the mag- 
netic flux density perpendicular to the two-dimensional 
lattice B = ph/(qe2a 2 ) is described by the number p/q 
of magnetic flux quanta li/e per plaquette 2a 2 . This 
differs from the random flux model studied previously 
i32ll . where the constant magnetic field part was absent. 
The random part is generated by the local fluxes <p XtZ , 
which are uniformly distributed -f/2 < <p x ~ < f/2 with 
zero mean. The disorder strength / can be varied within 
the interval from f/(h/e) = to f/(h/e) - 1. 



2.2. Orthogonal symmetry 

For the chiral orthogonal symmetry, we consider two 
different models. Both are defined on a square lattice. 
In the first, we set t z /to = 1 and take the disordered 
hopping terms as random real numbers that are defined 

as 



W 

t x = to exp — s, 

to 



(6) 



where {e} is a set of uncorrected random numbers with 
box probability distribution |e| < 1/2. The strength 
of the disorder was varied in the range W/to = 2 to 
W/to = 10. In the second model, the transfer terms in 
both the x and z direction are random numbers which are 
box-distributed about fo within the interval [to - W.,/2, 
to + W s /2] with possible disorder strengths < WJto < 
2. In the latter model, single sites may get isolated and 
decoupled from the remaining 2D lattice for W/to — > 2. 



2.3. Symplectic symmetry 

The off-diagonal disorder as given by Eq. (6) with 
W/to = 4 was used in the chiral symplectic model, 
where we consider a chiral version of the Ando model 
J5, 33]. The hopping terms t„„> are now 2x2 matrices 



t\\ - t 



t± = t x 



(7) 



where c 2 + s 2 = 1, s = 1/2, and the disorder in the 
hopping t x is given by Eq. ||6). 
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Figure 2: The bricklayer lattice (left) shares the topology of the hon- 
eycomb lattice (right). The red lines display the attached perfect semi- 
infinite leads used in the calculation of the scaling variables z, and the 
two-terminal conductance. 



2.4. Lattice topology 

The two-dimensional lattices of size L x L z on which 
the above models are defined, are either a regular square 
lattice or a bricklayer lattice which mimics the honey- 
comb lattice (Fig. |2]». The width L and the length L z of 
the lattice are taken to be similar L z L for the inves- 
tigations of the DOS and for the conductance, whereas 
quasi-one-dimensional samples with L z ~ 10 9 lattice 
units are used in the analysis of the scaling properties of 
the models. 

3. Numerical methods 

For the numerical calculation of the two-terminal 
conductance, we attach two ideal (without disorder) 
semi-infinite leads to the finite sample having width L 
and length L z and apply the well known Economou- 
Soukoulis formula 0411 
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where T is the corresponding transmission matrix and 
the X,- parameterize the eigenvalues of the hermitian ma- 
trix r"' T. The parameter A^ pen determines the number of 
open channels in the attached leads. Owing to the dis- 
order, the conductance is a statistical variable. There- 
fore, an ensemble of finite samples, which differ only 
in the microscopic realization of the disorder, is consid- 
ered and the mean value, the variance, and probability 
distribution of the conductance is evaluated using the 
algorithm described in 0~5ll . 

We define the dimensionless quantities z, = x,L/L- 
to be used for the scaling analysis in the vicinity of 
the critical point olSll . In the limit L z » L, the ratio 
Z\/{2L) converges to the smallest Lyapunov exponent 
y = lim£„^ a3 xi/(2L :; ) where x\ is the smallest positive 
eigenvalue of ln(T'T'). We calculated numerically the 
first two parameters z\ and zi for Q ID systems of length 
up to L z ~ 10 8 L in order to achieve a relative uncertainty 
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Figure 3: The spectrum of scaling variables \zt\ for the Q1D systems 
with chiral orthogonal symmetry. A lattice LxL. with square topology 
is considered with an odd value L = 65 a and two boundary conditions 
(BC) Dirichlet (left) and periodic (right) imposed in the transversal di- 
rection. For Dirichlet BC, the model is chiral orthogonal at the band 
center E = and zi(E = 0) = 0, in agreement with Eq. )10t . The 
symmetry changes to orthogonal for non-zero energy E. Note the de- 
generacy in the spectrum of scaling variables for the Dirichlet BC. 
This degeneracy is broken for non-zero energies E. There is no chiral 
symmetry for periodic BC. The disorder strength in the vertical hop- 
ping {4} is W = 5. Similar results for the chiral unitary symmetry 
were published in 12^1 . 



for z\ of order of 10 4 . The data obtained for z\{E,L) 
and zi(E, L) were fitted to scaling formulae (Q]) and (f2]). 

4. Properties of chiral systems 

In this section, we summarize our numerical results 
obtained for various disordered systems obeying chiral 
symmetry. 

4.1. Scaling variables Zi 

The calculated spectrum of scaling variables z, = 
XiL/L z is strongly influenced by the presence of the chi- 
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ral symmetry [20, 22]. They take on also negative values 
for chiral unitary systems. Our data J2stl confirmed that 



|Z2a-ll = \Z2a\ = C [fl - 1/2] 

for an even number of open channels, and 



\Zi\ = 



c Int[a/2] 
(c/2) [a- 1/2] 



Dirichlet BC 
periodic BC 



(9) 



(10) 



for an odd number of channels. From Eq. ( [TUt it follows 
that the mean value of the smallest zi is zero. 

The same relations hold also for chiral orthogonal 
systems. Figure [3] shows the spectrum of scaling vari- 
ables for iVopen odd and two types of boundary condi- 
tions. The numerical data indeed confirm that z\ =0 
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Figure 4: The smallest scaling variable zi for chiral unitary systems 
without a constant magnetic field, calculated on the bricklayer lattice. 
The number of open channels N p Cn is identical with the width of the 
layer. Periodic BC were used in the transversal direction. Note that 
the energy decreases up to E/to ~ 10~ 30 . Since quartic precision in 
the numerical calculations is necessary to treat such energies, we were 
able to calculate the scaling variables z, for these energies only for 
narrow samples (L < 24), due to the accuracy and CPU requirements. 



when Dirichlet BC are imposed. This is a reason for the 
non-Ohmic behavior of the mean conductance shown 
in Fig.Q] Also, Fig.[3]shows that the degeneracy of the 
spectrum of absolute values |z, | is lifted for any non-zero 
energy. A more detailed analysis of the scaling variables 
Zi for the chiral unitary system is given in Ref. 1E9I1 . 

4.2. Scaling analysis 

Scaling theory predicts that in the vicinity of the crit- 
ical point E — 0, the smallest scaling variable zi is a 
function of only one parameter, 



Zi(E,L)=F(x) 



(11) 



where x — InL/ ln(^(E)/^o), and £(£) is the correlation 
length 
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In 
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with unknown energy parameter £0 and exponent k. 

Figure [4] shows typical numerical data obtained in 
the vicinity of the chiral critical point. The data con- 
firm that both zi and zi can be described as a function 
of In E. However, a logarithmic energy dependence of 
£(£), given by Eq. ( fT2l can be observed only for very 
small values of the energy, typically |£7?ol < 10" 8 . 

For orthogonal and symplectic models as described 
in section |2] with disorder given by (O, we verified the 
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Figure 5: Scaling of the smallest scaling variable zi for chiral orthog- 
onal models with hopping terms t x given by Eq. {6). The disorder 
strength W/to is given in the Figure. The scaling variable x on the 
horizontal axis is given by Eqs. {1} and (2) with f o = 1, £o ~ ~4.5 to 
and the exponent k close to the value 2/3 (37h . The inset shows the 
disorder dependence of zi (E = 0) for L = 96 a. 



scaling behavior (flZb for each given disorder strength. 
In all systems we used periodic BC and even N open . 
This choice enables us also to compare systems with 
and without a constant magnetic field. Systems with 
odd N open and Dirichlet BC are numerically not accessi- 
ble since the smallest z, is zero for E — and very small 
for non-zero energy (Fig. [3}. 

We fit the numerical data for z\ and zi to the poly- 
nomial function of x an d extract the critical exponent 
k. We found that for all models k * 2/3, in a gree ment 
with previous theoretical predictions 1 24 , 3(| 31]. As 
an example, we show in Fig. [5] the result of our scal- 
ing analysis for the chiral orthogonal model and various 
strengths of the disorder. The disorder dependence of 
Z\(E = 0, W) is plotted in the inset. This shows that the 
the value of the disorder W represents an extra parame- 
ter that defines the respective model. 

Although the addition of a constant magnetic field 
to the chiral orthogonal model breaks the time-reversal 
symmetry, it does not influence the value of the expo- 
nent k as 2/3. The only exception is the chiral unitary 
model without constant magnetic field where k ^ 1/2 
i32ll (the raw data for this model are shown in Fig. |4|. 
More details about the scaling at chiral quantum critical 
points can be found elsewhere 13711 . 

We did not succeed in obtaining the logarithmic 
scaling from the raw data for the mean conductance 
J26L E9I1 . The reason is that for each particular sam- 



Figure 6: The probability distribution P(g - (g)) for the chiral uni- 
tary system without magnetic field and three different values of /, 
The width of the distribution varg as 0.187 does not depend on 
/. The solid line is the conductance distribution for the chiral uni- 
tary model with / = 0.5 h/e and with an additional magnetic field 
B = (l/l2)h/(ea 2 )(L/a = 72, mean conductance (g) = 2.82, varg = 
0.207). The inset shows the mean conductance (g) c for various val- 
ues of the random field /. Square samples of size 257 X 257 were 
considered with Dirichlet BC in the transversal direction. Note that 
the strength of the disorder is limited in the present model. Even for 
/ = 1 h/e, we get a rather large mean conductance, (g) = 1.45 e 2 //i. 



pie the conductance is given as a sum of contributions 
of all channels. Since the sum z\ + Zi is almost con- 
stant also for non-zero energies, the conductance de- 
pends only very weakly on the energy in the vicinity of 
the critical point. This tiny energy dependence, if any at 
all, cannot be extracted from our numerical data. 

4.3. Conductance 

At the critical point, the probability distribution of 
the conductance P(g) does not depend on the size of 
the system, and possesses the typical shape for a given 
universality class. In Fig. [6] we show the conductance 
distribution for the chiral unitary system without a con- 
stant magnetic field. The data confirm that the distribu- 
tion does not depend on the strength of the disorder /, 
which in this case measures the fluctuation of the flux 
through the plaquettes. Since the disorder is weak, the 
mean value of the conductance is large, and decreases 
when the disorder increases, as shown in the inset. The 
critical distribution is Gaussian with a universal width, 
varg = 0.187. 

Of particular interest is the comparison of the two 
unitary chiral systems that differ from each other only 
by the presence of a constant magnetic field. As dis- 



5 



1.5 



P(g) 



0.5 



p/q <g> var g 



• 





0.84 


0.18 


o 


1/16 


0.87 


0.23 


V 


1/8 


1.05 


0.25 



• V 

V 




-0.0001 

Energy E/t n 



0.0001 



0.0002 



Figure 7: The conductance distribution for the two-dimensional dis- 
ordered chiral orthogonal model on a square lattice with random hop- 
ping terms {6), and with a constant magnetic field, B = (p/q) h/ea 2 . 
The disorder strength is W/tg = 5, and the size of the system L/a = 64. 
The legend gives the ratio p/q which determines the magnetic field 
and corresponding mean conductance and the variance. Note the typ- 
ical non-analytical form of the distribution at g = 1 e 2 /h f38ll . 



Figure 8: The depression in the density of states p(E) near the quan- 
tum critical point at E = for the random magnetic flux model with 
(upper blue curve) and without (lower red curve) a spatial constant 
magnetic flux density defined on a bricklayer lattice. The histogram 
bin-width is about 4 X 10~ 6 . The dip in the DOS depends on both 
the disorder strength and the sample size. The brown and magenta 
solid lines are fitting curves l27l) . The RMF strength and system 
size are / = 0.5 h/e, L = L x = 96 a for B = 0, and / = 0.5 h/e, 
L = L x = 120 a in the case of B = (1/12) hj (lea 2 ), respectively. 



cussed in the previous section, these two systems ex- 
hibit a different critical exponent k. Surprisingly, the 
conductance distributions for these two models are al- 
most identical as shown in Fig. [6] Also, the mean con- 
ductance is larger than 1 even in the case of strongest 
disorder possible (/ = 1 h/e). Therefore, our chiral uni- 
tary model does not allow studies of strongly disordered 
chiral systems. 

Smaller values of the critical conductance can be ob- 
tained for systems with random hopping terms as given 
by Eq. ([6}. Figure|7]shows the critical conductance dis- 
tribution for system with real hopping terms given by 
Eq. (O with disorder W — 5 and various strengths of the 
constant magnetic field. While the system without mag- 
netic field belongs to the chiral orthogonal universality 
class, the constant magnetic field changes the symme- 
try to chiral unitary. In spite of the different symme- 
try classes, no significant difference between the criti- 
cal distributions is observed. Also the critical exponent 
k ss 2/3 is the same for both models 037I1 . 



in the DOS at the band center 112711 
.ST 
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4.4. Density of states 

In many 'standard' disordered models, the density of 
states can be considered to be constant in the critical 
region. A very strong energy dependence of the DOS 
can affect the accuracy of the numerical scaling. A typi- 
cal example is the metal-insulator transition in the band 
tails (see, for instance, J39]). Recently it was found 
that chiral unitary models can exhibit a very narrow dip 



It has been sug- 
that these 'microgaps' appearing in dis- 
ordered chiral models are the consequence of the non- 
perturbative ergodic regime. This can happen for times 
large compared with the diffusion time when the local- 
ization length exceeds the sample size. 

As an example, the DOS in the vicinity to the chi- 
ral unitary critical point is shown in Fig. [8] The en- 
ergy range of the DOS depression gets narrower with in- 
creasing RMF disorder strength for zero constant mag- 
netic field but it becomes broader in the case of a finite 



constant magnetic field 12711 . In the latter case, two ad- 
ditional quantum Hall critical points (QHCPs) showing 
the usual power-law divergence with critical exponent 
v =i 2.4 112611 exist symmetrically about the chiral point 
at E — 0. Within the achieved uncertainty this value 
for the chiral quantum Hall case obtained for a brick- 
layer lattice is compatible with the recent high-precision 
estimates v ^ 2.6 14 ll 14211 obtained for the Chalker- 
Coddington network model 114311 . With increasing dis- 
order, the QHCPs move further apart l26Tl and so they 
can be easily distinguished from the chiral one remain- 
ing at E = and studied here. 

To check whether this is a general feature of systems 
with chiral symmetry, we calculated numerically also 
the density of states for chiral orthogonal and for chiral 
symplectic models. The results are shown in Figs. [9] 
and [10] Similar to the case of chiral unitary symmetry, 
the chiral symplectic system shows a spectral gap in the 
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Figure 9: The density of states for the chiral symplectic model calcu- 
lated by diagonalizing 2 X 10 5 samples of size 16 X 16. A dip in the 
density of states is present only for weak disorder and small size of 
the lattice, but disappears when disorder and/or lattice size increases. 
The inset shows the density of states for disorder W = 2 and L = 1 6 
and L = 32. 
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Figure 10: The density of states for the chiral orthogonal model de- 
fined on a square lattice with real random hopping in both directions 
obtained by diagonalization of more than 2 X 10 4 samples of size 
80 X 80 and disorder strengths W s /t = 1.0 (blue) and WJt = 1.5 
(red). No dip but an extra peak is observed in the density of states 
close to the quantum critical point E = 0. 



DOS at E = 0. For a square system L — 64 and W — 3 
hopping disorder in both spatial directions, the narrow 
gap is still visible. With increasing disorder strength or 
system size, however, the gap becomes narrower until 
it seems to disappear due to the insufficient numerical 
resolution. 

In the chiral orthogonal case, instead of a dip, an extra 
peak is seen on top of the disorder broadened DOS at the 
quantum critical point E — 0. This result was obtained 
for a square lattice model with random hopping in both 
directions where van Hove singularities appear at the 
band center and not at E/to — ±1 as for the bricklayer 
(hexagonal) lattice. An explanation for the occurrence 
of the extra peak is the possible isolation of single lat- 
tice sites that become disconnected from the remaining 
system with increasing disorder. These sites contribute 
to the DOS with eigenvalues close to zero. 

5. Conclusion 

Two dimensional systems with chiral symmetry ex- 
hibit interesting new physical properties, not observed 
in standard disordered models. The most striking fea- 
ture is that the existence of a chiral quantum critical 
point is determined by the boundary conditions and by 
the parity of the size of the system. The chiral critical 
point coincides with narrow structures in the energy de- 
pendence of the density of states which depend both on 
the disorder strength and on the system size. To study 



the intriguing scaling behavior, a logarithmic scaling 
ansatz is necessary which replaces the usual power-law 
dependence applied in ordinary Anderson localization 
with diagonal disorder. 

Owing to the logarithmic energy dependence of the 
correlation length, the critical region is very narrow. 
Typically, it is narrower than the 10~ 8 th part of the band- 
width. This makes the scaling analysis of the chiral crit- 
ical behavior rather difficult. To achieve relevant data 
with sufficient accuracy, extremely long systems must 
be calculated, in some cases with quartic precision (up 
to 34 digits) of arithmetic operations. This is probably 
the reason why the logarithmic scaling was not observed 
in previous numerical studies |29[ 44-47]. Our numer- 
ical data indicate that all chiral systems can be divided 
into two classes: the first one consists of only one model 
with random magnetic flux and zero magnetic field. For 
this model, we observed that the critical exponent which 
governs logarithmic energy dependence of the correla- 
tion length is k 1 /2. For all other models, the critical 
exponent is k 2/3, independent of the symmetry. 

PM thanks Project VEGA 0633/09 for financial support. 
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